July 21, 1997 



MC-TH-97/12 



FROM DEEP INELASTIC SCATTERING TO 
PHOTOPRODUCTION: A UNIFIED APPROACH 

G. Kerley and G. Shaw 
Department of Physics and Astronomy, Schuster Laboratory 
University of Manchester, Manchester M13 9PL, U.K. 



ABSTRACT 

The strikingly different high energy behaviours of real photoabsorption cross-sections 
with = and the low x proton structure function at large are studied from 
a laboratory frame viewpoint, in which the x and dependence reflects the space- 
time structure of the interaction. This is done using a simple model which incor- 
porates hadron dominance, but attributes the striking enhancement observed at 
HERA at very low x and high to contributions from heavy long-lived fluctua- 
tions of the incoming photon. Earlier published predictions of the model for the 
then unknown behaviour of the structure function at small x and intermediate 
are shown to be strikingly confirmed by recent experimental data. A simultaneous 
analysis of real photoabsorption data and structure function data for < x < 0.1 
and < < 15 GeV^ is then reported. An excellent fit is obtained, with all 
parameters in the restricted ranges allowed by other physical requirements. 



1 Introduction 



At high energies, the real photoabsorption cross-section 

a^p{iy, Q2 = 0) ^ apz/"^~^ + ai?i/"«"^ 

is characterised by the intercept ap ~ 1.08 of the "soft" pomeron famihar from 
hadron scattering [Q, together with an effective contribution from lower lying Regge 
trajectories with ap ^ 0.55. The same behaviour characterises the proton structure 
function in the intermediate x-region 0.02 < x < 0.1, but experiments at HERA 
0, ^ in 1993 first observed the much sharper rise in the structure function at small 
X < 0.02 and large Q"^ > 8 GeV^ associated with the "hard" pomeron. Here we 
focus on the transition between these contrasting behaviours in the intermediate 
region < < 8 GeV^ using an approach suggested by one of us [Q. This was 
used to predict the behaviour at intermediate in 1994 but there was no data 
with which to confront the predictions. Subsequently the experimental situation has 
been transformed by precise data both at high and in the hitherto unexplored 
region of low x and intermediate Q^. |^, ^. Here we examine whether the new 
data at intermediate agree with the prior predictions of the model; and whether 
its greater precision and its extension to smaller x-values enable the intercept of the 
hard pomeron to be determined more precisely by the data. 

The model used is described in the next section. However it will be useful to first 
summarise the underlying approach This adopts the laboratory frame viewpoint 
typical of hadron dominance models, and emphasizes the role of the coherence length 

2t/ _ 1 1 1 

m2 + Q2 Mxl + myq^ ^ Mx ' ^ ^ 

which represents the typical distance travelled by a vacuum fluctuation of the photon 
of mass m. Our assumptions on the nature of the dominant states at different 
coherent lengths are summarised in Figure 1. At short coherence lengths, I < 
/o ~ 1 f, corresponding to large x values x > 0.1, we assume that these states are 
essentially bare qq pairs and a hadronic component of the photon has not developed. 
At moderate coherence lengths, they are assumed to have the single hadron-like 
behaviour expected for constituent qq pairs or vector mesons. However at very long 
coherence lengths, the photon fluctuations can eventually develop into the hadronic 
final states observed in e~^e~ annihilations^. For low masses, these states can be 
approximated by a sum of vector mesons; but for high masses complicated jet-like 
final states are observed. The idea is that the new phenomena at very small x are 
associated with complicated jet-like states which only play a role for masses m and 
coherence lengths I which are greater than some critical values m > rrij and / > Ic- 
By (|I]) this implies that they are confined to energies 

iy>iyc = m%/2 {Q^ = 0) (2) 



-'^For an attempt to understand the role of the coherence length from a more fundamental 
dynamical viewpoint, see Del Duca, Brodsky and Hoyer [pi. 
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in photoproduction, and to the kinematic region 

for structure functions. The remaining states with m < mj and/or / < Ic will 
completely dominate outside this region. They are included using conventional 
hadron dominance ideas which emphasize the link between real and virtual photons 
and provide a natural framework for describing the onset of scaling in the region^ 
xc < X < 0.1 . For example, they give a good description of the transition be- 
tween real photoabsorption and deep inelastic scattering data on protons in the 
intermediate x-region 0.02 < x < 0.1 0, 0; and they successfully predict the ob- 
served shadowing behaviour for real photoabsorption and deep inelastic scattering 
on nuclei in the same x region [[IT], |T^, |T^. 



2 A simple model 

In the hadron dominance mo deig, the total cross-section for photoabsorption is given 
by an expression of the form 



Q) = Jdm Jdm' ^ q,]^^,, ^ q,^ (4) 

corresponding to Figure 2. Conventional hadron physics and duality with the par- 
ton model both require large "diagonal" m = m! and large "off-diagonal" m ^ m' 
components to be present; and the latter must be explicitly incorporated if the ap- 
proximate scaling behaviour of the shadowing for deep inelastic scattering on heavy 
nucei is to be understood |Tl|, 0, |T3|. However, these off-diagonal contributions are 



dominated by terms in which the two mass values m and m' are not very different. 
Equation (^) is therefore often replaced by a simple diagonal approximation 

where p(z/, m?) is an effective quantity, meant to represent the more complicated 
structure on average^. Single hadron-like behaviour of the intermediate states can 
be incorporated by assuming a Regge-type energy dependence 

p^(j/,m2) = /p(m2)z/---i 



^In hadron dominance models which give approximate scaUng, the bulk of the contributions 
for a given comes from intermediate states with < rn^ < nQ^, with n a small integer. The 
average value is of order m? « Q^, for which the coherence length I « l/2Mx by (|^). 

•^For reviews, see for example Q and and references therein. 

^ If were exact, in the context of scaling in e+ annihilation and deep inelastic scattering 
it would imply that the total cross-sections am for states m scattering on nucleons decreased 
approximately as m~^. This counter- intuitive result is sometimes known as the Gribov paradox. 
It implies that the mean free paths in nuclear matter increase like m^, so that shadowing would 
die away at large at fixed x. These problems do not occur in "off-diagonal models" of the form 
when the effective quantity p{v,m?) results from destructive interference between diagonal 
and off-diagonal terms in (^.( Again, see |l^, |l^ and references therein.) 
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for the dominant diffractive contribution, where /p(m^) is a smoothly varying func- 
tion of mass chosen to lead to approximate scaling at large . In particular, if we 
assume 

(6) 
(7) 



this gives 



av 



dm 



(m 



2\l-Of 



(m2 + Q2)2 



for the soft pomeron contribution to the structure function, where 

Ma 



A 



This formula obviously embodies single hadron-like behaviour ap ^ 1.08 for all 
intermediate states. Here we modify it by adjusting the contributions from m > mj 
and I > Ici which is meant to roughly characterize the region in which single hadron- 
like behaviour has given way to a more complicated behaviour. This behaviour 
is unknown, so we parameterize it by the simplest possible generalization of (^), 
in which A, ap are replaced by new parameters A, a'p to allow for a different 
magnitude and energy dependence. In this way we arrive at a representation of the 
form 



+ Axu^p 



dm 



[m 



2\l-ap 



[m2 + Q2)2 

(m^V-^' 



[9{m'^j - m^) + e{m'^ - m'^j)e{lc - I)] 
e{m^ - m])e{l - Ic) . 



(8) 



For data fits, this is conveniently rewritten as 
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= 9{m\ 
+ Axv'^p 



m ; 



Axv 



dm 



dm^ 



[m 



2\l-ap 



[m2 + Q2)2 



fm2)i-"p 



+ e{m] - ml) 



m2 + Q2)2 

Axu'^^ 



+ Axu'^'' 



dm 



(m 



2\l-ap 



(m2 + Q2)2 



2^"''(^2 + g2)2 



(9) 



where the integrals can conveniently approximated by rapidly convergent series, as 
shown in the Appendix, and the critical coherence length is replaced by an effective 
critical mass, 

ml{x,Q')=Q'' ^ 



xMA 



- 1 



(10) 



which, unlike mj (and Ic), depends on Q"^ and x. The real photon equivalent is 



anv,{])=d{2v/lc-m]) 



ap — 1 



+ e{m] - 



m , 



Bu' 
{2u/l 



rrin 



im : 



-dp 



c) 



2u/l 



c) 



Bu 



ap~l. 



mn 
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where 

2.568Mpap' 2.568Mpa^' ^ ' 

and the numerical coefficients are chosen to give the cross section in mb for masses 
in GeV. 

In addition to the contributions associated with the hard and soft pomerons, 
there is also an additional contribution associated with lower lying Regge poles. 
This is a small correction in the region we are considering, and we incorporate it 
using the simple empirical form 

F^{x, Q') = Anx'--- (t^^) ' (13) 

VQ + Or/ 

proposed by Donnachie and Landshoff [|l5l, with the Regge intercept an = 0.55. 
The total structure function is then given by 

F2(a;,Q2) = Ff + (14) 

In the "pre-HERA" regions u < uc (Q^ = 0) and x > xciQ"^), (U) reduces to 
the simpler form (|^, and our parameterisation ( |T^ ) can be made numerically almost 
identical to the empirical Donnachie and Landshoff parameterisation |jl5|, which is 
known to give an excellent fit to photoproduction data and to structure function 
data in the intermediate region 0.1 > a; > 0.01. Correspondingly, data in this region 
effectively determines the parameters A,mo,AR and qr in (^ |13|). The additional 
parameters A, a'p, mj, Ic in (|^) describe the behaviour at very small x, with mj, Ic 
constrained to be a few GeV and a few fermi respectively. 



3 Predictions at intermediate 

In 1994 P|, the parameters A^ttiq^Ar and an were fixed by fitting the simple rep- 
resentation d^, |T3|, |n|)) to the "pre-HERA" data on real photoabsorption and the 
structure function at intermediate x values 0.02 < x < 0.1. The additional parame- 
ters A, a'p, mj, Ic in the final form (^ were then determined by extending the fit to 
include the then recently available HI data [It] at small x and large > SGeV"^; 



and used to predict the behaviour at intermediate {0 < < SGeV'^), where no 
data existed. In particular, the dramatic rise at small x observed at high was 
predicted to decrease rapidly and shift towards x = as becomes smaller. In this 
short section, we compare these predictions with the subsequent HI measurements 
0, in order to test the validity of the modelQ. 

Before doing this, we consider the quality of the 1994 fits for > 8 GeV^. This 
is illustrated in Figure 3 for the Regge intercepts 

ap = 0.55 ap = 1.08 , (15) 



^ Since the 1994 fits were based on the 1993 HI data at large Q^, we initially compare their 
predictions at lower to later data from the same experiment. The data from ZEUS and other 
experiments will be discussed explicitly in the next section. 
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together with both a "conventional" value a'p = 1.27 and a very unconventional 
choice a'p = 1.08 for the intercept of the long-hved, high mass states associated 
with the hard pomeron in this approach. As can be seen, this intercept was not well 
determined by the existing data, because of the limited x-range covered. However 
the 1996 HI data are more precise and extend to smaller x-values. Plotting 
the same curves against this new data in Figure 4 leads to a clear preference for 
a'p = 1.27, although there are small discrepancies between the precise new data and 
the fit, since the parameters of the latter were determined by data of lower accuracy. 

We now compare the 1994 predictions for the intermediate region < < 
SGeV^ with the 1996 HI data, restricting ourselves to the case a'p = 1.27 in the 
light of the previous discussion^. Since in 1994 there was no data at all in this region 
that could be used to help determine the parameters of the fit, the curves shown in 
Figure 5 represent a genuine prediction. When taken in conjunction with Figure 4, 
they show that the dependence predicted by the model is in excellent agreement 
with experiment. 



4 A global analysis 

In the last section we focussed upon the 1996 HI data at small x and intermediate 
because the parameters of the 1994 predictions were fixed by high Q"^ data 
from the 1994 HI experiment. Here we report the results of a simultaneous chi- 
squared fit to all the recent photoabsorption data in the wider kinematic region 
< a; < 0.1 andQ < < 15 GeV^. In doing so, we allow all the parameters in 
the diffractive term (|^), including the Regge intercepts ap and a'p, to be determined 
by photoabsorption alone, in order to see whether the resulting values accord with 



reasonable physical expectations. However, since the non-diffractive term (13) makes 
only a small contribution in this region, we keep the Regge intercept fixed at the 
value an = 0.55 obtained from hadron scattering data, leaving Aji and as variable 
parameters. 

The data used in the fit comprise the small x data of the HI 0, and ZEUS ^ 
collaborations; the small and intermediate x data of the E665 collaboration; the 



intermediate-x data of the NMC collaboration; and the real photoabsorption 



cross sections |T^, |T9|, |2T|, Obviously, we need to consider the consistency of 
the different experiments. The ZEUS and HI data cover a similar kinematic region, 
and are consistent within errors, as we shall see. In addition, for x ~ 0.01 there is 
data from both E665 and NMC at low and intermediate Q^, as well as high data 
from ZEUS and HI. As can be seen from Fig. 6, the experiments are consistent with 
each other, except for the three NMC points at (x, Q^) = (0.0125, 3.26) (0.0125, 4.52) 
and, especially, (0.008,3.47). 

^For 3 X 10~^ < X < 5 X 10"'^, the predictions are insensitive to the choice of a'p, but outside 
this range the value 1.27 rather than 1.08 is preferred. 

^We restrict ourselves to < 15 GeV^ since we are interested in the transition region from 
real photons to deep inelastic behaviour. Our simple model will of course break down eventually, 
since it gives exact scaling as ^ oo. 
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The resulting fits to the various data sets are shown in Figures 7-11 and the 
corresponding chi-squared contributions hsted in Table 1, where in all cases, the 
statistical and point systematic errors have been combined in quadrature. As can 
be seen, good fits are obtained for all data sets with the possible exception of the 
NMC data, which contributes 94 towards chi-squared from 66 data points. However 
it is clear from Fig. 10 that the fit is satisfactory, except for a few points, usually 
towards the edge of the kinematic range of the experiment, which make a very large 
contribution to . For example, the three NMC points discussed above contribute 
23 towards chi-squared; and on examining Fig. 6, it is difficult to see how any 
smooth curve, which fits the rest of the data, can do any better. 

With this small caveat, we conclude that the consistency of the data and the 
quality of the fit are very satisfactory. The corresponding parameter values are given 
in Table 2 together with their errors^ In addition, we have performed a series of 
seven fits in which one of the data sets listed in Table 1 is omitted, and the remaining 
six are fitted, in order check whether the parameter values are sensitive to small 
changes in the data set. The results are summarised in Table 3. As can be seen, the 
parameters are fairly stable against such changes, except for the parameters A^, 
associated with the small non- diffract ive term (0). These are only separately well 
determined^ if the NMC data, covering the intermediate x-region, are included. 

We now comment briefly on the values obtained for the parameters. In contrast 
to the 1994 flts, the intercept 

= 1.289 ±0.007 (16) 
of the hard pomeron is well determined, while the value 

ap = 1.059 ±0.007 

of the intercept of the soft pomeron is similar to, but perhaps slightly lower than, 
the value ap = 1.08 favoured by the analysis of hadronic total cross-sections [jl]] in 
terms of a simple Regge pole model of the pomeron. The value 

mo = 0.68 ±0.01 

is somewhat below the p mass, as predicted long ago by generalized vector dominance 
models with off-diagonal termsQ, while the value of the parameter 

mj = 2.51 ±0.02 

corresponding to the transition from resonance dominance to jet-like behaviour ac- 
cords well with e~ annihilation data. Finally the value for the characterictic 
distance over which jet-like behaviour develops for heavy states is 

Ic = 3.45 ± 0.06 / . 

^These errors are, of course, correlated; the full error matrix is available on request. 
^The ratio is well determined by real photoabsorption data. 

^"^This follows naturally from the destructive interference between the diagonal and off-diagonal 
terms (l4[ H). 



6 



This is also very reasonable, given the usual estimate / 1 f for the characteristic 
distance over which light qq states develop into vector mesons, based on, for example, 
the onset of nuclear shadowing effects in photoproduction on nuclei. 



5 Summary 

We have given a unified treatment of both real and virtual photoabsorption data in 
terms of a modified hadron dominance model, in which the striking enhancement 
observed at HERA at very low x and high is attributed to contributions from 
heavy long-lived fluctuations of the incoming photon. We have shown: 

• that the published predictions of the model for small x and low and interme- 
diate < 8 GeV^, where no data previously existed, are confirmed by recent 
data; 

• that the model gives an excellent fit to the much improved real photoabsorption 
and proton structure function data over the whole region < x < 0.1 and 
< Q2 < 15 GeV^] 

• that the values of the parameters of the model, determined from this real 
and virtual photoabsorption data alone, are in good agreement with physical 
expectations and other sources of information. 



6 Appendix: evaluation of integrals 

Data fits were performed using the MINUIT routine to minimise chi-squared with 
respect to the parameters of equation (H). This is computationally quite involved 
and speed and efficiency are at a premium. For this reason the integrals in (P) were 
evaluated by using a convenient series expansion rather than the normal methods 
of numerical quadrature. 

The integrals in (|^) can be easily written as linear combinations of integrals of 
the form 

^ dz ; < |e| < 1 (17) 



(Q2 + zf 

where the lower limits usually arise from the step functions, and can depend on x 
and Q^. Writing Q'^y = z, 



where 



and r = b/Q"^ 
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If r > 1, then 



J{r, e) = - ^— + - In 



If r < 1, then 



-e(l + e) V , \ , : (19) 
n(n + 1 + e)r("+i) J ^ ^ 



r^?/ = - / 7T^^?/- (20) 



Hence 



{1 + yY sin(7re) A/r (1 + ?/)^ 

^(r,e) = -^-J(l/r,-e) (21) 

and the series for r > 1 in (|T9|) can be used. These series have two advantages 
ensuring speed and accuracy of computation: 

• the series are alternating so, provided the number of terms is not too large in 
view of the numerical precision used, the error in the expansion is known and 
the degree of accuracy easily controllable; 

• the expansion variable is the ratio of the integrand lower limit to and for 
many of the points this is sufficiently far from unity for the series to converge 
adequately after only a few terms. Even for a unit ratio, the convergence is 
not prohibitively slow. 
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Tables 



Data source 


Number of points 




Structure Function data 


HI 


64 


35.7 




ZEUS 


47 


47.1 




NMC 


66 


93.9 




E665 


77 


74.0 


Real photo-absorption 


Caldweir73 


9 


8.8 




Caldwell'TS 


30 


33.8 




HERA 


3 


2.3 


All data 


296 


295.6 



Table 1: Breakdown of into contributions from different data sets to the global 



Diff ractive (pomeron) 



m5 
A 



m 



A 
Ic 



0.466 ± 0.008 
1.059 ± 0.007 
0.629 ± 0.022 

6.28 ± 0.08 
1.289 ±0.007 
0.457 ± 0.019 

17.5 ± 0.3 



Non- Diffr acti ve( Regge) 



Ar 

an 

OIR 



0.163 ± 0.029 
0.039 ± 0.009 
0.547(fixcd) 



Table 2: Parameter values and errors resulting from the global fit. (Natural units, 
based on the GeV , are used throughout.) 
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Omitted Data Set | 






Structure Function 


Photoabsorption 






HI 


/FITS 


F,66^ 
Ljyjyjo 


1 N iVl V. 


V^cLlLiVVcii / O 


V^cLlLLWcll 1 O 




2 

11 




465 


463 


508 


458 


463 


n /IRQ 


Up 








1 .UU ( 


1 rizLQ 


l.UOU 


1 064 

J- .UUt: 


A 


u.ozo 


0.628 


0.638 


0.655 


0.654 


0.653 


n 61 7 

U.Ol / 




5.823 


6.497 


6.303 


6.277 


5.945 


5.911 


6.339 


a'p 


1.291 


1.316 


1.290 


1.274 


1.286 


1.287 


1.290 


A 


0.432 


0.410 


0.454 


0.490 


0.458 


0.455 


0.451 


Ic 


13.48 


16.90 


17.42 


11.60 


17.82 


17.72 


17.32 


An, 


0.170 


0.164 


0.155 


0.012 


0.137 


0.138 


0.176 


<ih' 


0.0159 


0.0102 


0.0392 


0.00029 


0.0382 


0.0331 


0.0112 


x' 


259 


236 


221 


188 


286 


262 


293 


No. of Points 


232 


249 


219 


230 


287 


266 


293 



Table 3: Parameter and values for fits with one data set in turn excluded. 
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Figure Captions 



Figure 1 Hadronic behaviour of the photon for different values of the coherence 
length / and invariant mass squared of the photon fluctuations, where / = l/2Mx 
for m? = Q^. Here Ic and Mj are the critical values of the coherence length 
and fluctuation mass that separate the two types of hadronic behaviour. Another 
critical coherence length Zq ~ 1/ separates the hadronic and purely electomagnetic 
behaviour. 

Figure 2 The hadron dominance model Eq.(^). 

Figure 3 The 1994 fits § to the 1994 HI data ^ at = 12 GeV^, for the 
values a'p = 1.27(dashed line) and a'p = 1.08(sohd line). 

Figure 4 Comparison of the 1994 fits § with the 1996 HI data [§ at = 
12 GeV^. The curves again correspond to a'p = 1.27(dashed line) and a'p = 
1.08 (solid line). 

Figure 5 Comparison of the 1994 fits [|] with the 1996 HI data g] at = 1.5 
and 5 GeV^ for a'p = 1.27(dashed line). 

Figure 6 Comparison of the data from various experiments at x = 0.0125 (lower 
set) and x = 0.008(upper set), where the latter have been scaled by a factor of 
two, for clarity. The E665 points for x = 0.008 have been obtained by linearly 
interpolating between x = 0.007 and 0.009, and the HERA data have also been 
interpolated slightly from neighbouring points. The dashed(solid) lines show the 
result of the global fit described in the text for x = 0.0125(0.008) respectively. 

Figure 7 Comparison of the global fit with the 1996 HI data [^] for representative 

values. 

Figure 8 Comparison of the global fit with the 1996 ZEUS data 0. 
Figure 9 Comparison of the global fit with the E665 data ||^ for representative 
X- values. 

Figure 10 Comparison of the global fit with the NMC data Jl^ . 



Figure 11 Comparison of the global fit with the real photoabsorption data 
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Figure 1: Hadronic behaviour of the photon for different values of the coherence 
length / and invariant mass squared of the photon fluctuations, where / = l/2Mx 
for = Q^. Here Ic and Mj arc the critical values of the coherence length 
and fluctuation mass that separate the two types of hadronic behaviour. Another 
critical coherence length /q ~ 1/ separates the hadronic and purely electomagnetic 
behaviour. 



13 




14 




Figure 3: The 1994 fits § to the 1994 HI data |TB[ at = 12 GeV"^, for the values 
a'p = 1.27(dashed fine) and a'p = 1.08(sohd fine). 
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Figure 4: Comparison of the 1994 fits § with the 1996 HI data § at = 12 GeV'^. 
The curves again correspond to a'p = 1.27(dashed hne) and a'p = 1.08(sohd hue). 
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Figure 5: Comparison of the 1994 fits [§ with the 1996 HI data § at = 1.5 and 
5 Gel/2 ^ 1.27(dashed hne). 
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Figure 6: Comparison of the data from various experiments at a; = 0.0125(lower set) 
and X — 0.008 (upper set), where the latter have been scaled by a factor of two, for 
clarity. The E665 points for x — 0.008 have been obtained by linearly interpolating 
between x = 0.007 and 0.009, and the HERA data have also been interpolated 
slightly from neighbouring points. The dashed(solid) lines show the result of the 
global fit described in the text for x = 0.0125(0.008) respectively 
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Figure 7: Comparison of the global fit with the 1996 HI data for representative 
values. 
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Figure 8: Comparison of the global fit with the 1996 ZEUS data 0. 
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Figure 9: Comparison of the global fit with the E665 data 
X- values. 
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Figure 10: Comparison of the global fit with the NMC data . 
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